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^^ Abstract. A recent paper of Arnold, Falk, and Winther [Bull. Amer. Math. 

O I Soc. 47 (2010), 281-354] showed that a large class of mixed finite element 

^ methods can be formulated naturally on Hilbert complexes, where using a 

^^ Galerkin-likc approach, one solves a variational problem on a finite-dimensional 

subcomplcx. In a seemingly unrelated research direction, Dziuk [Lecture Notes 
in Math., vol. 1357 (1988), 142—155] analyzed a class of nodal finite elements 
" for the Laplace— Beltrami equation on smooth 2-surfaces approximated by a 

^_^ piecewise-linear triangulation; Demlow later extended this analysis [SIAM 

_^ J. Numer. Anal, 47 (2009), 805-827] to 3-surfaces, as well as to higher- 

/^ order surface approximation. In this article, we bring these lines of research 

\^ together, first developing a framework for the analysis of variational crimes 

• in abstract Hilbert complexes, and then applying this abstract framework to 

i -~H the setting of finite element exterior calculus on hypersurfaces. Our framework 

jrt extends the work of Arnold, Falk, and Winther to problems that violate their 

asubcomplex assumption, allowing for the extension of finite element exterior 
calculus to approximate domains, most notably the Hodge-de Rham complex 
on approximate manifolds. As an application of the latter, we recover Dziuk's 
^». ^ and Demlow's a priori estimates for 2- and 3-surfaces, demonstrating that 

K^ surface finite element methods can be analyzed completely within this abstract 

framework. Moreover, our results generalize these earlier estimates dramatically, 
extending them from nodal finite elements for Laplace-Beltrami to mixed finite 
elements for the Hodge Laplacian, and from 2- and 3-dimensional hypersurfaces 
to those of arbitrary dimension. By developing this analytical framework 
using a combination of general tools from differential geometry and functional 
l/j analysis, we are led to a more geometric analysis of surface finite element 

C_,^ methods, whereby the main results become more transparent. 
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1. Introduction 

The aim of this paper is to bring together three distinct ideas that have influenced, 
in separate ways, the development and analysis of geometric finite element methods 
for elliptic partial differential equations. 

The first idea is that of a variational crime. Suppose we have a variational 
problem of the form: Find u £ V such that 

(1) B{u,v) = F{v), VweF, 

where F is a Hilbert space, B: V x V ^ M. is a, bounded, coercive bilinear form, 
and F G V* is a bounded linear functional, li Vh C V is a subspace (usually 
finite-dimensional), then one can obtain an approximate solution by solving the 
Galerkin variational problem: Find u^ G Vh such that 

B{uh,v)^F{v), \fveVh. 

This is the typical abstract setting for finite element methods. However, for many 
problems of interest, especially finite element methods on surfaces or on domains 
with curved boundaries, one cannot efficiently compute the bilinear form B (•, •) or 
the functional F (•) on a subspace of V. Instead, one must take an approximating 
space Vh 'f. V, along with an approximate bilinear form Bh : V/i x V/i -> M and 
functional Fh G V^*, and formulate the generalized Galerkin variational problem: 
Find Uh € Vh such that 

(2) Bh{uh,v) = Fh{v), yveVh. 

Such modifications to the original variational problem are called "variational crimes." 
There is a well-understood framework for the analysis of a large class of variational 



GEOMETRIC VARIATIONAL CRIMES 3 

crimes, represented by the Strang lemmas [7j: for instance, the first and second 
Strang lemmas allow for the complete analysis of numerical quadrature, the use 
of geometric modeling technology such as isoparametric elements, and many other 
examples of variational crimes. 

The emergence of surface finite elements represents a second distinct idea that 
has influenced the development of geometric finite element methods. The analysis 
of surface finite element methods, which by construction are "criminal" methods, 
has required a more sophisticated approach that exploits the specific nature of the 
crime in order to obtain a satisfactory error analysis; this custom-tailored analysis 
contrasts with the more general approach given by the Strang lemmas. The surface 
finite element research area was effectively initiated with the 1988 article of Dziuk 
]W\ , although there is related work appearing about ten years earlier by Nedelec [IS] . 
While there was some activity in the area during the 1990s (cf. |T71[TT]), beginning 
in 2001 there was a tremendous expansion of research in the general area of surface 
finite element methods, with many applications arising in material science, biology, 
and astrophysics; examples include [M l [T0 l [T2 l [T3 l [T9 l [T8 l fT5 t fT4 ] . 

The third distinct idea that has had a major influence on the development of 
geometric methods is that of mixed finite elements, whose early success in areas such 
as computational electromagnetics was later found to have surprising connections 
with the calculus of exterior differential forms, including de Rham cohomology and 
Hodge theory [6l [27l [28l [2^ . This has culminated, very recently, in the powerful 
theory of finite element exterior calculus developed by Arnold, Falk, and Winther 
[H [3] . A key insight of the latter work, from a functional-analytic point of view, is 
that a mixed variational problem can be posed on a Hilbert complex: a differential 
complex of Hilbert spaces, in the sense of Briining and Lesch [5]. Galerkin-type mixed 
methods are then obtained by solving the variational problem on a finite-dimensional 
subcomplex. 

In this article, we bring these lines of research together, first developing a 
framework for the analysis of variational crimes in abstract Hilbert complexes, 
and then applying this abstract framework to the setting of finite element exterior 
calculus on hypersurfaces. Our framework extends the work of Arnold, Falk, and 
Winther [3] to problems that violate their subcomplex assumption, allowing for 
the extension of finite element exterior calculus to approximate domains, most 
notably the Hodge-de Rham complex on approximate manifolds. As an application 
of the latter, we recover Dziuk's [TB] and Demlow's [T3] a priori estimates for 2- 
and 3-surfaces, demonstrating that surface finite element methods can be analyzed 
completely within this abstract framework. Moreover, our results generalize these 
earlier estimates dramatically, extending them from nodal finite elements for Laplace- 
Beltrami to mixed finite elements for the Hodge Laplacian, and from 2- and 3- 
dimensional hypersurfaces to those of arbitrary dimension. By developing this 
analytical framework using a combination of general tools from differential geometry 
and functional analysis, we are led to a more geometric analysis of surface finite 
element methods, whereby the main results become more transparent. 

The remainder of the article is organized as follows. In [Section "2} we review 
the abstract framework of Hilbert complexes, which plays a central role in the 
work of Arnold, Falk, and Winther [31 on finite element exterior calculus. This 
includes a brief introduction to Hilbert complexes and their morphisms, domain 
complexes, Hodge decomposition, the Poincare inequality, the Hodge Laplacian, 
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mixed variational problems, and approximation using Hilbert subcomplexes. In 
[Section 3} we consider the approximation of a Hilbert complex by a second complex, 
related to the first complex through an injective morphism rather than through 
subcomplex inclusion. Since this morphism is not necessarily unitary (i.e., inner- 
product preserving), this allows the approximating complex to have a different 
inner product, which only approximates that of the original complex. We develop 
some basic results for the pair of complexes and the maps between them, and then 
prove error estimates for generalized Galerkin-type approximations of solutions to 
variational problems using the approximating complex; these estimates generalize the 
results of Arnold, Falk, and Winther !3j to "external" approximations. Our results 
may be viewed as establishing Strang-type lemmas for approximating variational 
problems in Hilbert complexes. Finally, in |Section~4{ we apply the framework 
developed in [Section 3] to the Hodge-de Rham complex of differential forms on a 
compact, oriented Riemannian manifold. We first review Hodge-de Rham theory, 
and then consider a pair of Riemannian manifolds related by diffeomorphisms, 
establishing estimates for the maps needed to apply the generalized Hilbert complex 
approximation framework. After reviewing the concept of a tubular neighborhood, 
we then consider the specific case of Euclidean hypersurfaces. We subsequently 
show how the results of the previous sections recover the analysis framework and a 
priori estimates of Dziuk |16) , Demlow and Dziuk 15J, Dcmlow 14 , and moreover 
extend their results from scalar functions on 2- and 3-surfaces to general fc-forms on 
arbitrary dimensional hypersurfaces. We also indicate how our results generalize 
the a priori estimates of Dziuk |16) , Demlow |14j from nodal finite element methods 
for the Laplace-Beltrami operator to mixed finite element methods for the Hodge 
Laplacian. 

2. Review of Hilbert complexes 

In this section, we quickly review the abstract framework of Hilbert complexes, 
which forms the heart of the analysis in Arnold, Falk, and Winther [3] for mixed 
finite element methods. Just as the space of L^ functions is a prototypical example 
of a Hilbert space, the prototypical example of a Hilbert complex to keep in mind is 
the i^-de Rham complex of differential forms. (This example will be discussed at 



greater length in Section 4 ) After stating the basic definitions, we will summarize 
some of the key results from Arnold, Falk, and Winther [3] on mixed variational 
problems and their numerical approximation using Hilbert subcomplexes. The 
interested reader may also refer to Briining and Lesch [8^ for a comprehensive 
treatment of Hilbert complexes from the viewpoint of functional analysis. 

2.1. Basic definitions. Let us introduce the basic objects of study, Hilbert com- 
plexes, and their morphisms. 

Definition 2.1. A Hilbert complex {W,d) consists of a sequence of Hilbert spaces 
W'', along with closed, densely-defined linear maps d*^ : V'' C W''' — !• V'^^^ C W'^^^, 
possibly unbounded, such that d*^ o d*^"^ = for each k. 

• • • ), yk'l > yk ), yk+l s. . . . 



This Hilbert complex is said to be bounded if d*^ is a bounded linear map from W'' 
to W''^^ for each k, i.e., (M^, d) is a cochain complex in the category of Hilbert 
spaces. It is said to be closed if the image d^V^ is closed in M^'^+i for each k. 
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Definition 2.2. Given two Hilbert complexes, (M^, d) and {W',d'), a morphism 
of Hilbert complexes f : W ^>- W consists of a sequence of bounded linear maps 
fk-iYk^ ^r'k guch that f'^V'' C V"" and d"'/'' = f^+^d'' for each k. That is, the 
following diagram commutes: 

. . . ). yk !, yfc+1 > ■ ■ ■ 



■ ■ ■ > V'^ > yik+i ). . . . 

By analogy with cochain complexes, it is possible to define notions of cocycles, 
coboundaries, harmonic forms, and cohomology spaces for Hilbert complexes. 

Definition 2.3. Given a Hilbert complex (VF, d), the space of k-cocycles is the 
kernel 3^^ = kerd*"', the space of k- coboundaries is the image 'B*'' = d''^^V^~^, 
the kth harmonic space is the intersection Sj'' — '^'^ Ci 05*^''""', and the kth reduced 
cohomology space is the quotient 3'^/'8'^. When *B'^ is closed, 3'^/23'^ is simply caUed 
the kth cohomology space, and is identical to reduced cohomology. 

Remark 1. One can show that the harmonic space io*^ is isomorphic to the reduced 
cohomology space 3'^/^'^- For a closed complex, this is identical to the usual 
cohomology space 3'^/S'^, since S*^ is closed for each k. 

Definition 2.4. Given a morphism of Hilbert complexes f : W ^ W, the induced 
map on (reduced) cohomology is defined by [z] i— > [fz] , where [z] denotes the (reduced) 
cohomology class of the cocycle z. 

In general, the differentials d of a Hilbert complex may be unbounded linear 
maps. However, given an arbitrary Hilbert complex (W^, d), it is always possible to 
construct a bounded complex having the same domains and maps, as follows. 

Definition 2.5. Given a Hilbert complex {W, d), the domain complex {V, d) consists 
of the domains V'^ C W'^ , endowed with the graph inner product 

{u,v)yk = {u,v)^yk + {d''u,d''v)^yk+i- 

Remark 2. Since d'^ is a closed map, each V'^ is closed with respect to the norm 
induced by the graph inner product. Also, each map d'' is bounded, since 

Thus, the domain complex is a bounded Hilbert complex; moreover, it is a closed 
complex if and only if {W, d) is closed. 

For the remainder of the paper, we will follow the simplified notation used by 
Arnold, Falk, and Winther [3]: the VF-inner product and norm will be written simply 
as (•,•) and ||-|j, without subscripts, while the y- inner product and norm will be 
written explicitly as {■,-)y and ||-||y. 

2.2. Hodge decomposition and Poincare inequality. The Helmholtz decom- 
position states that a rapidly-decaying vector field on K^ can be decomposed into 
curl- free and divergence- free components, i.e., the vector field can be written as the 
sum of the gradient of a scalar potential and the curl of a vector potential. For 
differential forms, this is generahzed by the Hodge decomposition, which states that 
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any differential form can be written as a sum of exact, coexact, and harmonic com- 
ponents. Here, we recall an even further generalization of the Hodge decomposition 
to arbitrary Hilbert complexes; this immediately gives rise to an abstract version of 
the Poincare inequality, which will be crucial to much of the later analysis. 

Following Briining and Lesch [8], we can decompose each space W'' in terms of 
orthogonal subspaces, 

where the final expression is known as the weak Hodge decomposition. For the 
domain complex {V,d), the spaces 3*^, *B'', and io'"' are the same as for (PF, d), and 
consequently we get the decomposition 

T/fe = S^©io'=©3'=-Lv^ 

where 3'^"'"^ = 3'^"'""' HV''. In particular, if {W, d) is a closed Hilbert complex, then 
the image QS'^ is a closed subspace, so we have the strong Hodge decomposition 

W^ = «B'=®i3''©3''^"'', 
and likewise for the domain complex, 

From here on, following the notation of Arnold, Falk, and Winther [3], we will 
simply write 3^^^ in place of 3'^^*' when there can be no confusion. 

Lemma 2.6 (Abstract Poincare Inequality). If {V,d) is a bounded, closed Hilbert 
complex, then there exists a constant cp such that 

||«||y<cp||d'=z;||^, Vve^"^. 

Proof. The map d'^ is a bounded bijection from 3*^^ to *B*^+^, which are both 
closed subspaces, so the result follows immediately by applying the bounded inverse 
theorem. D 

Corollary 2.7. If (V,d) is the domain complex of a closed Hilbert complex (W,d), 
then 

IMIv<cp||dS||, Vi;e3'^. 

We close this subsection by defining the dual complex of a Hilbert complex, and 
recalling how the Hodge decomposition can be interpreted in terms of this complex. 

Definition 2.8. Given a Hilbert complex {W, d), the dual complex {W* ,d*) consists 
of the spaces Wj^ = W'^, and adjoint operators d^. = (d*""^)* : V^ C W^ -J> ^k-i ^ 

■ ■ ■ < v,u ^^^ V,* 4^ v,\, < 

Remark 3. Since the arrows in the dual complex point in the opposite direction, 
this is a Hilbert chain complex rather than a cochain complex. (The chain property 
d^ o d^_|_2 = follows easily from the cochain property d'"' o d*^"^ — 0.) Accordingly, 
we can define the k-cycles 3^ = kerd^ = sgfe-Lw and k-boundaries W^ = '^l+i^k ■ 

2*. 



The kth harmonic space can then be rewritten as Sj'' — 3'° H 3^; we also have 
■^k _ fg*J-w^ and thus 3'^^^''' — 05^. Therefore, the weak Hodge decomposition can 
be written as 



GEOMETRIC VARIATIONAL CRIMES 7 

and in particular, for a closed Hilbert complex, the strong Hodge decomposition 
now becomes 

2.3. The abstract Hodge Laplacian and mixed variational problem. To 

obtain a "mixed version" of the familiar Poisson equation —Am = / for scalar 
functions, we now follow Arnold, Falk, and Winther [3] in defining an abstract 
version of the Hodge Laplacian for Hilbert complexes. The abstract Hodge Laplacian 
is the operator L = dd* + d*d, which is an unbounded operator W'' — ?> W^ with 
domain 

DL^{ue V^ n T4* I du £ V7+1, d*u e v''-^} . 
If u G Dl solves Lu = f, then it satisfies the variational principle 

{du,dv) + {d*u,d*v) = {f,v), Vwe V'^nV^. 

However, as noted by Arnold, Falk, and Winther [3], there are some difficulties in 
using this variational principle for a finite element approximation. First, it may be 
difficult to construct finite elements for the space V^ H VjT . A second concern is the 
well-posedness of the problem. If we take any harmonic test function v £ S)^ , then 
the left-hand side vanishes, so {f^v) = 0; hence, a solution only exists if / _L j^'^. 
Furthermore, for any 9 G io'^ = 3'^ H 3^, we have dg = and d*g = 0; therefore, if u 
is a solution, then so is w + q. 

To avoid these existence and uniqueness issues, one can define instead the following 
mixed variational problem: Find {(T,u,p) E V'^^^ xV^ x S)'^ satisfying 



(3) (da,t;) 



Here, the first equation implies that a — d*u, which weakly enforces the condition 
u £ V'' f] V^ . Next, the second equation incorporates the additional term (p, w), 
which allows for solutions to exist even when f ^ S^^ . Finally, the third equation fixes 
the issue of non-uniqueness by requiring u ^- S)^ . The following result establishes 
the well-posedness of the problem ([S]). 

Theorem 2.9 (Arnold, Falk, and Winther [5], Theorem 3.1). Let {W,d) he a closed 
Hilbert complex with domain complex {V, d) . The mixed formulation of the abstract 
Hodge Laplacian is well-posed. That is, for any f G W^ , there exists a unique 
{a,u,p) € V^~^ xV^ X S^'' satisfying (Is]). Moreover, 

lklly + hlly + l|p||<c||/||, 
where c is a constant depending only on the Poincare constant cp in \Lemma 2.6\ 

To prove this, one observes that ^ can be rewritten as a standard variational 
problem (fTl) on the space V''~^ x V'^ x ^'^ , with the bilinear form 

B {a, u, p;T,v,q) = {a, r) - (u, dr) + {da, v) + {du, dv) + {p, v) — {u, q) 

and functional F(T,v,q) = {f,v). The well-posedness then follows immediately 
from the following theorem, which establishes the inf-sup condition for the bilinear 
form B. 



{a, t) - {u, dr) = 0, 


W e ¥>" 


{du, dv) + {p, v) = (/, 1;) , 


VveV'' 


{u,q)^Q, 


Vq e ro^ 
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Theorem 2.10 (Arnold, Falk, and Winther [3;, Theorem 3.2). Let {W,d) be a 
closed Hilhert complex with domain com,plcx (V,d). There exists a constant 7 > 0, 
depending only on the constant cp in the Poincare inequality (Lemma 2.0^, such 



{(^h,r) ~ {uh,dT) = 


= 0, 


Vr e Vt 


{duh,,dv) + {ph,v) = 


= (/,«), 


yv e v/:, 


{uh,q) = 


= 0, 


Vg e iot 



that for any {a, u,p) £V'' ^ xV'' x 9)'^, there exists (r, v, q) G V' ^ xV'^ xS^^ with 
B {(T,u,p-T,v,q) > 7 (||a||^, + ||u||^, + ||p||) (||t||^ + ||t;||^ + \\q\\) . 

From the weh-posedness result, it follows that there exists a bounded solution 
operator K : W^ — >■ W^ defined by Kf — u. 

2.4. Approximation by a subcomplex. In order to obtain approximate numeri- 
cal solutions to the mixed variational problem ([3]), Arnold, Falk, and Winther [3] 
suppose that one is given a (finite-dimensional) subcomplex Vu C V of the domain 
complex: that is, Vj^ C V^ is a Hilbert subspace for each fc, and the inclusion 
mapping ih'. Vu ^^ V is a, morphism of Hilbert complexes. By analogy with the 
Galerkin method, one can then consider the mixed variational problem on the 
subcomplex: Find {(Jh,Uh-,Ph) S Vl^~^ x V,^ x Sj'j^ satisfying 



(4) {dah,v) 



For the error analysis of this method, one more crucial assumption must be 
made: that there exists some Hilbert complex "projection" nh- V —>■ Vh- We put 
"projection" in quotes because this need not be the actual orthogonal projection 
i^ with respect to the inner product; indeed, that projection is not generally a 
morphism of Hilbert complexes, since it may not commute with the differentials. 
However, the map tt^ is V-bounded, surjective, and idempotent. It follows, then, 
that although it does not satisfy the optimality property of the true projection, it 
does still satisfy a quasi- optimality property, since 

\\u-TThu\\y^ inf \\{I ~ TTh) {U - V)\\y < \\I - TThW luf \\u ~ v\\y , 

veVh veVh 

where the first step follows from the idempotence of tt/i, i.e., {I — TTh)v = for 
all V a Vh- With this framework in place, the following error estimate can be 
established. 

Theorem 2.11 (Arnold, Falk, and Winther [3], Theorem 3.9). Let (T4,d) he a 
family of subcomplexes of the domain complex {V, d) of a closed Hilbert complex, 
parametrized by h and admitting uniformly V -bounded cochain projections, and let 
{a, u,p) e T/'^-i xV^xS)'' be the solution of ^ and {ah, Uh,Ph) S V/^''^ x V^ x S)l 
the solution of problem (HI). Then 

\W ~^h\\v + \W-Uh\\v + \\v-Ph\\ 

<C{ inf ||cr-r|U+ inf ||u — u||y+ inf ||p— g|L + M inf II-Pqs" ^ ^'ll v)i 
revl^-^ vev^; qev^ vev^ 



where ^ ~ ^^ ~ sup (/ — tt^) r 



||r||=l 



Therefore, if Vh is pointwise approximating, in the sense that vaiy^y^ \\u — w|| — > 
as /i — >■ for every u &V , then the numerical solution converges to the exact solution. 
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3. Analysis of variational crimes 



In this section, we extend the results of Arnold, Falk, and Winther [3], summarized 
in the previous section, by removing the requirement for Vh to be a subcomplex 
of V. The key point of departure is in the map ih'- Vh ^^ V; rather than being an 
inclusion, we require only that it is an injective morphism of Hilbert complexes, with 
the property that tt/,, o i^ is the identity. (The latter requirement simply corresponds 
to the earlier condition that tt^ be idempotent in the case of subcomplexes.) After 
stating some basic results about complexes equipped with such maps, we develop 
error estimates for the mixed variational problem and eigenvalue problem on Vh- 
These estimates contain two additional error terms, in addition to those in the 
analysis of Arnold, Falk, and Winther ^. These extra terms, analogous to those 
in the Strang lemmas for generalized Galerkin methods, measure the "severity" of 
two variational crimes: first, how well the right-hand side ihfh approximates /; and 
second, the extent to which i^ fails to be unitary. 

3.1. Approximation by an arbitrary complex. In order to approximate a 
Hilbert complex (VK, d), suppose we have another Hilbert complex (T/F/i,d/j), along 
with a pair of morphisms: an injection i^ : Wh ^-> W and a projecti on tt^ : W -^ W h, 
such that TT^ o ijj is the identity on Wj^ for each k. Recall that, by Definition 2.2 of 
a Hilbert complex morphism, the maps ijj and tt^ must be bounded for each k. The 
relationships among the domains and maps are illustrated in the following diagram: 



fe+i 



-^y= >v 



>v,^^^v/:^' >■■■. 

Arnold, Falk, and Winther ^3 consider the case where Wh C M^ is a subcomplex, 
and ih is the inclusion of Wh into W. In this special case, ih is unitary (i.e., an 
isometry), since for all u,v G Wj^, we have {ihU,ihv) = {u,v) = {u,v)i^. Indeed, 
if ih is unitary, then we can simply identify Wh with the subcomplex ihWh C W. 
More generally, though, we will consider cases where Wh (f. W, and where ih is not 
necessarily unitary. 

We begin by demonstrating some basic facts about these approximations. 

Theorem 3.1. If {W,d) is a bounded Hilbert complex, then so is (Wh,dh)- 

Prnnf lld'^ll — IItt'^+^j'^+^H'^ II — IItt'^+^H'^?'^ II < IItt'^+HI lld*^ II II?''' II -^ on H 

rrooj. ||Q,j|| — ||7r^ z^ a^|| — ||7r^ a i,^\\ ^ ||7r^ ||||a \\\\ih\\ < oo. u 

Theorem 3.2. If {W,d) is a closed Hilbert complex, then so is {Wh,dh). 

Proof. Assume that {W, d) is closed, so that each coboundary space S'^ is closed 
in W''. Now, since ih is a morphism, if Vh € Sf^ then ihVh € S*^, so *B^ C z^^*S'^. 
Conversely, since iTh is a morphism, if ihVh G S'' then Vh — i^hihVh G S|, so 
— isgfc ^ sgfc_ Therefore, *Bjj = i^^S'^, and since ih is bounded (and hence 
continuous), it follows that *B^ is closed. D 

Since tt^ o *^ = id^*:, this composition induces the identity map on the reduced 

cohomology space 3{5/*Bj^; thus ih induces an injection on reduced cohomology, 
while TTh induces a surjection. We now show that, given a certain approximation 
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condition on the harmonic spaces Sj , these induced maps are in fact isomorphisms 
(which are inverses of one another, since their composition is the identity) . 

Theorem 3.3. Let {W,d) and {Wh,dh) be Hilhert complexes, with morphisms 
hi ■ VVfi ^-J- W and tt/^: W —>■ W^ such that tt^ o *^ = id^^/fc for each k. If, for all k, 

||g-i^7r^g|| < \\q\\ , VqE^'', q ^ 0, 
then TTfi (and thus ih) induces an isomorphism on the reduced cohomology spaces. 

Proof. Since iTh induces a surjection on reduced cohomology, it suffices to show 
that this is also an injection. That is, given z G 3*^ with nhZ £ 03^, we must 
demonstrate that z S QS'^. Using the weak Hodge decomposition, write z = q + b, 
where q G Sj'^ and b E ^B'^. By assumption, tt^z G *B^, and since tt/,, is a morphism, 
TTh.b G *B{J as well. Thus, nhq = i^hZ — Hhb G S^, and since ih is also a morphism, 
ihT^hq E S*^ -L -^'^- Therefore, iuT^hq -L 9, which implies that q violates the inequality 
above, so we must have q — Q and hence z G 55*^. D 

Corollary 3.4. // [W, d) and {Wh, d/j) are closed Hilbert complexes, with morphisms 
TTfi and ih satisfying the above assumptions, then tt/j (and thus ih) induces an 
isomorphism on cohomology. 

Remark 4. This result is slightly more general than Arnold, Falk, and Winther [3l 
Theorem 3.4], which only treated the case of a bounded, closed Hilbert complex. 
However, the proof is essentially identical. 

Next, suppose that (V^,d) and {Vh,dh) are bounded, closed Hilbert complexes; 
for example, they may be the domain complexes corresponding, respectively, to 
closed complexes (W, d) and {Wh,dh). We now show that the Poincare inequality 
for (yh,dh) can be written entirely in terms of the Poincare constant for (V,d), 
denoted by cp, along with the operator norms of ih and tt/j. 

Theorem 3.5. Let {V,d) and (Vh,dh) be bounded, closed Hilbert complexes, with 
morphisms ih'-Vh^^V and tt/j : V ^ Vh such that tt^ o i^ = idyt for each k. Then 

Proof. Given Vh E 3|^, let z E 3^^^ be the unique element such that dz — dih.Vh — 
ihdhVh- Then, applying the abstract Poincare inequality on V, 

\\z\\v < Cp \\<iz\\y = Cp \\ihdhVh\\v < cp||i|;.^^|| ||d,jU/,.||y^ . 

It now suffices to show ||f/i||y < \Wh\\ ll^lly- Observe that Vh — tt/^z G Vj^, and 
furthermore, 

dhTThZ = TT/idz = TThihdhVh = dhVh, 

so Vh — TThZ G 3/5 -L Vh- Therefore, 

\\vh\fv^ = {vh,TThz)y^ + {vh,Vh-nhz)y^ = {vh,nhz)y^ < \\vh\\y^ \\nf^\\ \\z\\y , 

and the result follows. D 

Corollary 3.6. //(V,d) and (Vhjdh) are the domain complexes corresponding, 
respectively, to closed Hilbert complexes (W,d) and (Wh,dh), then 

\\vh\\y, < cp||41|||«n WdhVhWh^ ^Vh E 3^. 
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Finally, given the importance of the projection morphism tt^ in finite element 
exterior calculus, we now prove a short but useful result on the existence of such 
projections. In particular, the next theorem states how a projection morphism on 
another complex, W , can be "pulled back" to obtain one on W, as pictured in the 
following diagram: 




In [Section 4[ this will allow us to obtain a projection morphism for the de Rham 
complex on a manifold, by pulling back the usual projection defined on its piecewise- 
linear triangulation. 

Theorem 3.7. Let (VF,d) and (W/i,d/i) he Hilbert complexes with an injection 
morphism i^: Wh ^-> W. Suppose there exists another complex (M^',d') and a 
morphism f : W ^ W , such that i'f^ = f o ih' Wh ^^ W' is injective and has a 
corresponding projection morphism tt^ : W' — >■ Wh with ttJ^ o i'^ — idvi/^ • Then there 
also exists a projection morphism tt/i : M^ — >■ Wh such that nh o ih — idiy^ • 

Proof. Take nh = Tr'h o /• Then iTh o ih ^ n'f^, o f o ih ^ Ki ° */i = idv^^- □ 

3.2. Modified inner product and Hodge decomposition. As noted in the 
previous section, this generalized framework introduces some new complications, 
due to the possible non-unitarity of ih- The following result shows that the subspace 
ihWh C W can be identified with Wh, endowed with a modified inner product 
( J/i', •)^ instead of (•, ■)i^. This defines a modified Hilbert complex, which will be 
denoted by {ilW,dh)- 

Tiieorem 3.8. Let ih : Wh ^^ W be a morphism of Hilbert complexes, and define 
Jfc = if^*ik -w^^W^ for each k. Then 

{JhUh,Vh)h = {ihUh,ihVh) , yuh,Vh e Wh, 
is an inner product, which defines a Hilbert space structure on Wj^ . 

Proof {ihUh,ihVh) = {ih^hUh,Vh) h = {JhUh,Vh)h- This is an inner product, since 
ih is linear and injective. Moreover, Wj^ is closed with respect to the induced norm, 
since |jJhf/i|| < \\ih\\ W^hW and ih is bounded, so this is indeed a Hilbert space. D 

Remark 5. We use the notation Jh due to the similarity with the Jacobian determi- 
nant used in the "change of variables" formula for integration. Note that, although 
each JJl : Wj^ — > Wj^ is a bounded linear map, Jh is not necessarily a Hilbert complex 
automorphism. This is because, in general, d commutes with ih but not with its 
adjoint i^^. Also, clearly ?{j is unitary if and only if JJl = id^yfe. 

Now, if ih does not preserve the inner product, in particular it does not preserve 
orthogonality: that is, Uh -L Vh does not imply ihUh -L ihVh- This has significant 
implications for the Hodge decomposition, since although W}^ — *8^ ® i5^ ® 3^"'""''' 
is PlOi-orthogonal, it is generally not ijj H^-orthogonal. Therefore, we define the new, 
modified subspaces 

sj'h' = {ze 3f. I ihz ^ ^hK} > 3^"" = {« e Wl: I thv ± ihi'h} • 
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This gives a modified Hodge decomposition Wj^ = S^ © Sj'j^ © 3^^'^ , which is no 
longer necessarily W/i-orthogonal, but is now zJ^W^-orthogonal. As before, this also 
gives a modified Hodge decomposition for the domain complex V^ = *8^©io5f ©3/5^'- 

3.3. Stability and convergence of the mixed method. Let (VF, d) be a closed 
Hilbert complex with domain complex (V,d). To approximate a solution to the 
mixed variational problem ^, suppose that {Wh,dh) is another Hilbert complex 
with domain complex {Vh,dh), and that we have morphisms ih- V^ ^-)- V and 
TT/i : V" — >■ Vh such that T^h°'''h — i'^v'" ^^^ ^^*^h ^- We assume that i^ is W- 
bounded, so that it also can be extended to W^ ^-> W, but that tt/j might only be 
^-bounded. Then consider the solution of the following mixed variational problem: 
Find {ah,Uh,ph) e V^~^ x V^ x S)l satisfying 

(ct/i, Th)f^ - {uh, dhVi)^ ^ 0, Vt,i e ^/f"\ 

(5) {dh(Jh,Vh)h + {<^hUh,dhVh)h + {Ph,Vh)h = {fh,Vh)h : V-y/i G V^ , 

This corresponds to the generalized variational problem (pi) with bilinear form 

Bhicrh,Uh,Ph;Th,Vh,qh) = {o-h,Th)h- {uhAhTh)h 

+ {diiah^Vh),^ + {dhUh,dhVh)h + {ph,Vh)h - {'^h,qh)h 

and functional F^ {Th,Vh,qh) = {fh,Vh)h- ^^^ following theorem establishes the 
inf-sup condition for the mixed method ([5]). 

Theorem 3.9. Let (V,d) be the domain complex of a closed Hilbert complex (W,d), 
and let (V/j, d/i) be a family of domain complexes of closed Hilbert complexes (W^, d/j), 
equipped with uniformly W -bounded inclusion morphisms i^ : Vh ^~> V and V -bounded 
projection morphisms tt/i : V —^ Vh satisfying nf^ o i^ = idyk . Then there exists a 
constant jh > 0, depending only on cp and the norms of ih and tt^, such that for 
any {ah, Uh,Ph) & V/^~^ x V^ x ij^, there exists (r/,, Vh,qh) e V^~^ x V^ x Sj^ where 

Bh io-h,Uh,Ph;Th,Vh,qh) 

> Ih (IkftJIv^ + \\uh\\v^ + \\Ph\\h) {hhWv^ + \\vh\\v^ + hhWh) ■ 



Proof. This is just Theorem 2.10 applied to the Hilbert complex {Vh, dh), combined 



with the fact that the Poincare constant is cp Wi^hW ||i/i|| by Theorem 3.5 D 



Remark 6. Since we have assumed that the morphisms ih and tt/j are uniformly 
bounded with respect to h, it follows that the inf-sup constants 7/1 can be bounded 
below by some constant, which is independent of h. 

The goal, for the remainder of this section, will be to control the error 

Ik - ihC^hWv + II" - ihUhWv + \\P - ^hPhW , 
where {a,u,p) is a solution to pi) and {ah,Uh,Ph) is a solution to ([5]). To do this, 
it will be helpful to introduce the following modified mixed problem on ihV' Find 
K,<,p'J e Vt' X V/: X ^'^^ satisfying 

{JhCTh, Th)h - {Jhu'h, dhTh)h = 0, Vr,, e Vl^~^, 

(6) {Jhdha'h,Vh)h + {Jhdhu'h,dhVh)h + {Jhp'h,Vh)h = {ihfiVh)h , Vw,, G V^ , 



{Jhulq'h)h=0, Vg;, eio 



/k 
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This has the corresponding bihncar form 

+ {JhAh(T'h,Vh)f^ + {Jh<^hu'i„dhVh)h + {Jhp'h,Vh)h - {Jhu'h,q'h)h > 

and the functional F;^ {tu, vh, q'^) = {ilf, w/i),,- 

This is precisely equivalent to the mixed problem on the subcomplex ihVh C V, 
which has the bounded cochain projection ih o tth: V ^f ihYh- Therefore, the 
stability and convergence analysis of Arnold, Falk, and Winther p] can be applied 
immediately to this modified discrete problem. In the end, we will obtain the desired 
bound by applying the triangle inequality, 

(7) \\cT -ih(7h\\v + \W -ihUhWv + \\P- ihPhW 

< \\(T - ihcr'hWv + ll" - ^hu'hWv + lb - ihP'hW 

+ \\ih{(yh-(yh)\\v + \\ih{uh~u'^)\\v + ¥h{Ph -p'h)\\ ■ 
Observe that, since ih is bounded, we can write 

\\ih {<yh - <y'h)\\v + 11*'^ K - Oily + \\ih {Ph~p'h)\\ 

< C{\\(Th -<y'h\\vy, + \Wh~u'h\\v^ +\\Ph-p'h\\h) ' 

so it will suffice to control the error between solutions to ^ and ^ in V^. 

Theorem 3.10. Under the assumptions of Theorem 3.S\ suppose that {o'h,Uh,Ph) G 
V^^^ X y,f X ioj; is a solution to ^ and {a'^, u'hiP'h) ^ ^h^^ ^ ^h ^ -^/f *-' "' solution 
to ^. Then 

\Wh - 't'hWv, + h'^ - <\\v, + \\Ph-Ph\\k < c iWh - tUWh + \\i- JhW 11/11) • 

Proof. For any {T,v,q) G V/^^^ x V/^ x i^^, we can write 

Bh{<^h-T,Uh-v,Ph -q]Th,Vh,qh) = B^ {(Jh - f^'h^Uh - u'f^.ph - p'h]Th,Vh,qh) 

+ Bh Wh - T, M/i ~v,p'i^- q; Th, Vh, qh) ■ 
Ignoring the first term momentarily, observe for the second term that 

Bh{(^h^ u'h , Ph ; Vi , Vh ,qh) = B'f-^ {(t'^ , u'h , p'h \ Th , Wft , qh ) 

+ ((/ - Jh) (T'h,rh)h - {{I - Jh) u'^AhTh)h + {{I - Jh) Aha'y^,Vh)h 

+ {{I - Jh) <^hu'h, dhVh)h + {{J - Jh)Ph,Vh)h - {{I - Jh) u'h,qh)h , 
so by the variational principles ^ and ([6|, 

B'h (o-/i, u',^,p'h',Th,Vh, qh) = {i*hf, Vh)h - {Jhu'h,qh)h , 
Bh{(Th,Uh,Ph\Th,Vh,qh) = {fh,Vh)h- 
Therefore, 

Bh (o-/i - o-^, Uh - u'f^,Ph - p'h, Th, Vh, qh) = {fh ~ i*hf, Vh)h + {u'h,1h)h 

- {{I - Jh)<7'h,Th)h + {{I - Jh)u'h,dhTh)^- {{I- Jh)dhcr'^,vh),^ 

- {{I - Jh)dhu'f^AhVh)h- {{I - Jh)p'h,Vh)h, 
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SO using the boundedness of the bilinear form and Cauchy-Schwarz, we get the 
upper bound 

Bh{(^h-T,Uh-v,Ph -q-:Th,Vh,qh) 

< c{\\h - *UII„ + WPs^u'hWh + \\i- All (IkJJIv, + IKIIv., + Wp'hh) 

+ hh-'^Wv,, + h'h~v\\y^ + \\p',^ -'7llh)(lk'.|ly;. + \\vh\\v,, + ll'?'i||/i)- 
Next, [Theorem 33] gives the lower bound 

BhicTh -T,Uh- v,ph -q;Th,Vh,qh) 

> Ih {\Wh ~ T\\y^ + \\uh - v\\y^ + \\ph - q\\,^) {\\Th\\v^ + \\vh\\v^ + ||g/J|,,) 

for some {Th,Vh,qh) G l^/f^^xV'^f x-Q^' where 7/1 can be bounded below independently 
of h. Therefore, combining the upper and lower bounds and dividing out ||T;i||y + 
hhWv^ + \\qh\\h, we get 

\Wh - T\\y^ + \\uh - v\\y^ + \\ph - q\\f^ 

< C{\\h - ^U\\,, + \\P^,yhh + \\I- JhW (Ikklly, + IKIIv^. + WPhh) 

+ h'h ~ T\\y^ + IK - t;||,,^^ + Up',, - g||,,J . 

This expression can be simplified considerably by choosing r = crjj, v — u'j^, and 
q = Ps^i^p'fi, so applying the triangle inequality gives the error estimate 

hh-(T'h\\v^ + \\uh-u'^\\y^ + \\ph -p'hWh 

< c {\\h - ^UWh + ll^i-,.<llh + 11^ - M 11/11 + Wp'h - q\\h) ■ 

All that remains is to deal with the terms H^fl^w'^jH,, and \\p'f^ — q\\y^. First, since u'f^ is 
j^T^-orthogonal to i^'^f , the modified Hodge decomposition lets us write u'f^ = u'^+m^, 
where w'^g G S^ and m^ G 3j5^'- Now, observe that Psj^u't^ = since *8j^ _L i^^, and 
furthermore P^^Jhu'j^ ~ since w^ € 3^^"""' implies Jtu'^ J- 3^- Therefore, 

ll^-^.<lk = ll^io.^ll. - ll^>5. (^ - AXII, < C 11/ - J,|| 11/11 . 

Next, since p'^^ e S)'l^ C 3^, the Hodge decomposition gives p'f^ — P<ShP'h + P-<Oh.P'h = 
P^hP'h + 9- Also, similar to the previous term, since p'f^ G Si^ we have JhP'h -L '^\^ 
so P<s,^JhP'h = 0- Thus, 

Wp'h ~q\\h^ WP'sJhWh = II^'b. U - A)pkL < ^ 11^ - All ll/ll . 

Therefore, these two terms can be combined with the existing ||/ — J,j|| ||/|| term, 
leaving the final error estimate, 

hh - o'^v, + II"'. - <K + \\ph-p'h\\H < G {\\h - *UII„ + 11^ - All ll/ll) , 

as desired, which completes the proof. D 
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Corollary 3.11. If{a,u,p) e V''^^ xV'^ xSj'' is a solution to ^ and {ah,Uh,Ph) G 
V^'^^ xV/^x iojj is a solution to ^, then 

Ik - ihO-hWv + II" - ihUhWv + \\P - ^hPhW 

<C{ inf Ik — T||y+ inf ||m — w||t/+ i^f \\P^ qWv+I^ ^^^ \\P<su - v\\y 
reihVl^-^ veihV^ qeihV^ veihV^ 

+ \\h-^yh + \\I-Jh\\\\f\\), 

where fi = fif^ = sup (/ — i^^Tr^^) r . 

\\r\\=l 

Proof. Use the triangle inequality, as in ([7]), and then apply Theorem 2.11 and 



Theorem 3.10 to bound the respective error terms. D 



This theorem establishes convergence, as long as our approximations satisfy 
||/ — JfiW — > and ||/;i — ihfWh ~^ ^ when h -^ 0. This raises the question of how to 
choose fh E V,f ; although clearly fh — i*yj will work, in many cases this cannot be 
computed efficiently. The next result demonstrates that, if II/j : W^ —J- W^ is any 
bounded linear projection (i.e., satisfying II/j o i^ = lA-^k), then simply choosing 
fh = ^hf is sufficient to get a quasi-optimally convergent solution. 

Theorem 3.12. // II/j : W'' —>■ Wj^ is a family of linear projections, bounded 
uniformly with respect to h, then we have the inequality 

||^,J-^;J||,<C(||/-A||||/||+ inf \\f-^). 
Proof. Using the triangle inequality, we write 

ll(n, - ii) /ii, < WiUf, - iit^^u^,) /ii, + IK*;, - ^l^^,n^,) f\\,^ 

= \\iI-^l^h)Il,J\\,^ + \\^l{I-^hU,,)f\\^ 

<\\i-Jh\\\\ii,,.fh + \m\\ii-z,,n„)f\\ 

<C(||/-A||||/||+ inf 11/ -011), 

where the final step follows from the M^-boundedness of Hh and the quasi-optimality 
property of J - z,,n,,, i.e., (/ - ^H^) f = {I - ih^h) (/ - 0) for any <j) e ihWJ^. D 

3.4. Remarks on obtaining improved error estimates. Arnold, Falk, and 
Winther [3] were also able to obtain improved error estimates by making some 
additional assumptions: namely, that tt/j is VF-bounded rather than merely V- 
bounded, and that the Hilbert complex V satisfies a certain compactness property. 
With these assumptions, the continuous solution operator K: W^ — >■ W^ becomes 
a compact operator, and hence converts the pointwise convergence of / — tt/j — > 
(which follows from the quasi-optimality property) to norm convergence. This norm 
convergence is essential for applying the so-called "Aubin-Nitsche trick" (also known 
as "i^ lifting"), where one obtains improved estimates by applying the solution 
operator to the error term itself. Roughly speaking, one needs norm convergence, 
rather than pointwise convergence, since the solution operator is being applied to 
quantities that depend on the parameter h. 

However, there are no such improved estimates for the additional error terms 
obtained in the previous subsection. Essentially, this is because norm convergence 
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is already required for ||/ — J/j|| — > as /i — > 0, and there is no analogous quasi- 
optimality result for J/j as there is for vr/j. Therefore, these terms remain the same, 
and the improved estimates only apply to the terms already analyzed by Arnold, 
Falk, and Winther [3 for the subcomplex case. 

3.5. Convergence of the eigenvalue problem. While we have primarily focused 
on the numerical approximation of the mixed variational problem, Arnold, Falk, and 
Winther [3] also analyzed an eigenvalue problem associated to the Hodge Laplacian. 
The extension of their eigenvalue convergence result to non-subcomplexes is fairly 
straightforward, and follows from the results already given in this section, as we 
will now show. 

Consider the eigenvalue problem 

((T, r) - (u, dr) = 0, Vrey''^-\ 

(8) {da, v) + (du, dv) + {p, w) = A (m, u) , Vu G V ^ 

{u,q)=0, Vqei3^ 

the discrete problem 

{<^h,Th)h - {uh,dhTh)h = 0; ^'^h e V,f"\ 

(9) {dhCTh,Vh)h + {dhUh,dhVh)h + {Ph,Vh)h^ \h{uh,Vh)h, Vu/i G V^'', 

and the modified discrete problem 

( Acr^: '^hjh - {Jhu'h, dhTh),^ = 0, Vr/i e V^''"\ 

(10) {Jhdh(T'h,Vh)i^ + {Jhdhu'i^AhVh)h + {Jhp'h-,Vh)h = K {'^'h^'"h)h : Vw/i G Vj^, 

(A<,<z;)^ = 0, yq'H^^h- 

As shown by Arnold, Falk, and Winther [3l Theorem 3.19], solutions to the subcom- 



plex problem (10) converge to those of (|8|, which follows immediately from the fact 
that ijiK'f^Pii converges to K in the C (W^ , W''^ operator norm. We now show that 
this result also holds for the problem Q. 

Theorem 3.13. Let {V,d) he the domain complex of a closed Hilbert com,plex {W,d) 
satisfying the compactness property, and let {Vh,dh) be a family of domain complexes 
of closed Hilbert complexes {Wh,dh), equipped with morphisms iu'- Wu ^^ W and 
TTh'- W ^ Wh such that tt^ oi^ — idyyk, where ih and tt/i are hounded uniformly with 
respect to h. Then the discrete eigenvalue problems Q converge to the problem ([8]). 

Proof. It suffices to show that ikK^Ph converges to K in the £ (W*^, W'^) operator 
norm. (As stated by Arnold, Falk, and Winther [^ , the sufficiency of norm conver- 
gence follows from Both, Brezzi, and Gastaldi [51.) Using the triangle inequality, we 
write 

\\K - ihKhPhW <\\K~ ihKPhW + Wh (K ~ Kh) PhW . 
The first term on the right-hand side converges to zero, by Arnold, Falk, and 
Winther [2 Corollary 3.17]. For the second term, recall that ih and tt/j are assumed 
to be bounded uniformly with respect to /i, and since \\Ph\\ — \\''^hPihWh\\ ^ IKft-ll; 
it follows that P^ is bounded uniformly with respect to h, as well. Therefore, it 
suffices to control ] j iiT^ — iiT/jJI in C{W^,WJl). However, the earlier analysis in 



GEOMETRIC VARIATIONAL CRIMES 17 



Theorem 3.10 shows that ||A'^ — Kh\\ < C\\I — Jh\\, which completes the proof of 



convergence. D 

4. Application to differential forms on Riemannian manifolds 



In this section, we apply the framework developed in Section 3 to the Hodge- 
de Rham complex of differential forms on a compact oriented Riemannian manifold. 
We will begin by first recalling the basic definitions of the de Rham complex of 
smooth forms; its completion as a Hilbert complex, called the L^-de Rham complex; 
and the corresponding domain complex, which dovetails with the theory of Sobolev 
spaces. Next, we discuss the general problem of approximating the de Rham 
complex on a manifold M by a family of "nearby" manifolds M^, each equipped 
with an orientation-preserving diffeomorphism tph : Mh — > M . We subsequently 
establish the correspondence between this setup and the generalized Hilbert complex 
approximation framework of |Section 3[ obtaining estimates for the appropriate maps, 
as needed. We then specialize the discussion a bit further by considering the case 
when M is a submanifold of some larger manifold N] in this case, the approximating 
submanifolds Mh C N can be taken to lie in a tubular neighborhood of M, and 
Lph '■ Mfi — )> M is obtained by projection along normals. 

Finally, we then look at the specific case where N = M" , and where we wish to 
approximate a solution on some to- dimensional Euclidean hypersurface M C M", 
n = 771 + 1, by finite elements defined on a piecewise- linear mesh M^ C M". This is 
now the realm of of surface finite element methods, as analyzed in Dziuk |16j . Demlow 
and Dziuk ^15j , Demlow I14j . We subsequently show how our results of the previous 
sections recover the analysis framework and a priori estimates of Dziuk [16] , Demlow 
and Dziuk [TS], Demlow [T3], extending their results from scalar functions on 2- 
and 3-surfaces to general fc-forms on arbitrary dimensional hypersurfaces. We also 
indicate how our results generalize the a priori estimates of Dziuk |16j . Demlow 
[l4] from nodal finite element methods for the Laplace-Beltrami operator to mixed 
finite element methods for the Hodge Laplacian. 

4.1. A brief reviev^r of Hodge— de Rham theory. Given a smooth, 777-dimensional 
manifold M, let H'^{M) denote the space of smooth fc-forms on M for /c = 0, 1, . . . , m, 
and let d'' : ^^{M) -> f7''+^(M ) be the exterior derivative for /c = 0, 1, . . . , 777 - 1. 
Then (r2(M),d) is a cochain complex, 

— > n°{M) -^ n\M) -^ S i7"(Af) — > 0. 

called the de Rham complex on Af . 

Suppose that, in addition, M is oriented and compact, and has a Riemannian 
metric g. Then, we can define the LP' -inner product of any u,v ^ n'^{M) to be 



("'«)l20(M) = / UA*gV^ {{U,V)) fig. 

^ ' Jm Jm 

Here, -kg : i}''{M) — )■ ri™~'^(Af ) is the Hodge star operator associated to the metric, 
((•, •)) : i}''{M) X rt''{M) -^ C°°{M) is the pointwise inner product induced by the 
metric and /ig is the Riemannian volume form. (The Hodge star is defined precisely 
so that uA-kgV = {{u,v)) Hg, and it follows that *g is an isometry.) The Hilbert 
space L^ri'"'(Af) is then defined, for each k, to be the completion of J7'^(Af) with 
respect to the L^-inner product. 
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To show that this forms a Hilbert complex (L^r2(Af), d), we must now define tlie 
weak exterior derivative d'" on some dense domain oi L^n^{M). Given u e L'^n'^{M), 
we say that w £ _L^ri''"+^(A/) is the weak exterior derivative of u, and write du = w, 
if 

where il^+^(Af) denotes the space of smooth {k + l)-forms with compact support. 
Therefore, one defines the dense domains H^^{M) C L'^VL^{M), consisting of 
elements in L?^^{M) that have a weak exterior derivative in L^O''+^(M). Thus, 
we have 

> Hn°{M) -^ H9}{M) -^ S HVl"^{M) > 0. 

where each HVt^{M) can be given the graph inner product 

("'«)m(M) = ("'«)L2n(Af) + (dW'dt;)^2ij(M) ■ 

(Note the similarity with the definition of the Sobolev spaces i/^, iJ(curl), and 
H (div).) Since each HQ''{AI) is complete, it follows that d'' is a closed operator; 
therefore, (L^ri(M),d) is indeed a Hilbert complex, and {Hn{M),d) is the corre- 
sponding domain complex. Furthermore, it can be shown that (L^n(M),d) satisfies 
the compactness condition, so these Hilbert complexes are in fact closed and satisfy 
the conditions necessary for the improved error estimates. (For more details on the 
construction of these complexes, see Arnold, Falk, and Winther [5].) 

4.2. Diffeomorphic Riemannian manifolds. Let {M,g) be an oriented, com- 
pact, rn-dimensional Riemannian manifold, and suppose {Mh,gh) is a family of 
oriented, compact Riemannian manifolds, parametrized by h and equipped with 
orientation-preserving diffeomorphisms Lph '■ Mh — > M . Now, since the pullback 
if*j^ : ^{M) — T' O {Mh) and pushforward Lp^^ : O (Af/j) — > VL{M) commute with the 
exterior derivative, they give a cochain isomorphism between the smooth de Rham 
complexes fi {Mh) and VL {M). 

We now show that these maps are bounded, and hence can be extended to Hilbert 
complex isomorphisms between LP'VL{Mh) and L'^Q,{M)^ following the results of 
Stern |29]. Given any point x G Mh, let {ei, . . . , Cm} be a positively-oriented, gh- 
orthonormal basis of the tangent space T^Mh, and let {/i, . . . , /„} be a positively- 
oriented, ^-orthonormal basis of T^^(a;)M. Then, with respect to these bases, the 
tangent map T^^h'- T^Mh — > T^p^^(^^)M can be represented by an to x to matrix 
$. Moreover, since Lpt is a diffeomorphism, the matrix $ has to strictly positive 
singular values, 

Oii{x) > ■ • ■ > am{x) > 0. 

These singular values are orthogonally invariant, so they are independent of the 
choice of basis at each x and iph{x). Hence, they are an intrinsic property of the 
diffeomorphism, and thus we refer to them as the singular values of iph at x. 

Theorem 4.1 (Stern |29], GoroUary 6). Let {Mh,gh) and {M,g) be oriented, 
vn- dimensional Riemannian manifolds, and let (ph '■ Mh — > M be an orientation- 
preserving diffeomorphism with singular values ai{x) > • • • > an{x) > at each 
x e Mh- Given p,q G [l,oo] such that 1/p + 1/q — 1, and some k = 0, . . . ,to., 
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suppose that the product (ai • • • am-k) {oim-k+i • ■ ■ Q^m) is bounded uniformly 

on Mh. Then, for any uj G L^Vt^ (Mh), 

II ^ \^/<] ( N^l/PII^-'- II I 

ll(ai---afe) '^(afe+i---a,„) '^||^ \\uj\\^ 

< \\iph*uj\\p < II (ai • ••a„_fc)^'^(a„_fc+i ■ ■ ■ a^y^ \ 



oo\M\p- 



Sketch of proof . At each point, a /c-form is /c-linear and totally antisymmetric. 
Therefore, the puUback is controlled pointwise by the product of the k largest 
singular values of (ph, while the pushforward is controlled by the product of the k 
largest singular values of tp^^ (i.e., the reciprocals of the k smallest singular values 
oi (fih)- Thus, we obtain pointwise inequalities 

\v*hV\ < Oil ■ ■ ■ ak i\ri\ ° ^h) , \^h*^\ < [{am-k+i ■ ■ ■ (^.^y |a;|] o (^^^ 

For the L^ upper bound, we can apply the pushforward inequality to get a factor 
of {a,n-k+i ■ ■ ■ c(„i) ^ in the integrand. Using the change of variables theorem 
introduces the Jacobian determinant ai • • • a™, so multiplying by this gives a factor 
of «! • • • am~k {oim-k+i ' ' ' ctm) ^ ■ We Can then use Holder's inequality to pull 
out the L°°-norm of this expression, and raising to the exponent 1/p gives 

II c \^/p ( \~1+1/PII 

||(ai • ■ -am-fc) [Om-k + l ■ ■ ■ OL,n) ||^ 

II I \1/P / \^l/<7 II 

= ||(ai •• -am-fe) "^(a„_fe+i •• •«,„) ' ||^, 

as desired. The lower bound follows in a similar fashion, starting with the identity 
w = (p*yLph*'^ and applying the pointwise pullback inequality. D 

Since M and M^ are compact, the uniform boundedness hypothesis of this 
theorem is clearly satisfied. Therefore, taking p = q = 2, it follows that the diffeo- 
morphism ^ph induces Hilbert complex isomorphisms ph* ■ L^^ (Mh) — > L^n{M) 
and ifl : L^n (M) ^ L'^Vt (il4). 

Now, take W — LF'Q,{M), and suppose that we have discrete subcomplexes Wh C 
L?'^{Mh) with inclusion morphisms i'^ : Wh ^-> L'^^{Mh), as well as projection 
morphisms ttJ^ : L^fl {Mh) — >■ Wh bounded uniformly in h. Following the approach 
of |Theorem 3.7[ we can pull these back to obtain the injection morphisms ih = 
Vh* ° i'h ■ ^h ^^ W and projection morphisms tt/j — n'l^ °V*h'- W ^ Wh, which satisfy 
T^h^ih = idwfc • An important consequence of this is stated in the following corollary 
of [Theorem 3.71 and [Theorem 4.11 

Corollary 4.2. Orientation-preserving diffeomorphisms induce an equivalence of 
families of finite element subcomplexes of the L^-de Rham complex with bounded 
cochain projections. In particular, any triangulation Th ^ M gives corresponding 
V^ andVr families (cf. Arnold, Falk, and Winther \^\3^) of piecewise-polynomial 
differential forms on M . 

Finally, let us see how this definition of ih can be used to control the error term 



/ — Jh\\- Theorem 4.1 implies that, for any Vh € V^, we have the estimate 



1/ \i/2 / ^-l/2||-l II I 

\[ai---ak)' (ak+i---am) \\^\\'"h\\h 



< \\ihVh\\ < II (ai •• -am-k)^ {am^k+i ■ ■■am) H^o ll^'^IL ' 
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and since Jh — i*hih, this implies 
\\ai ■ ■ ■ ak [ak+i ■ ■ ■ am)~ ||^ ll'^/ilU 



< \\JhVh\\h < ||"1 • • -Oim-k (ttm-fc+l • • -Oim) |L \\vh\\h 



This bounds the spectrum of the self-adjoint operator J/i, so finally we obtain a 
bound on the error term ||/ — J^W in terms of the singular values, 

ll-^- JhW < maxi -■ II - - ' - - ^"iii 

(11) ' 



1 - ||ai • • -afe (afe+i • • ■ a^n) 

1 - llai • • • am~k {a-m-k+i ■ ■ ■ am) 



-ii 



It follows that, if each singular value satisfies |1 — a^j < Ch"^^ , then ||/ — J/j|j < 
Ch'^'^^ as well, and moreover this will hold for every fc = 0, . . . , to. Obtaining such 
bounds on the singular values, for particular choices of Lpjn will be the topic of the 
next subsection. 

4.3. Tubular neighborhoods and EucHdean hyper surfaces. Suppose that 
(iV, 7) is an oriented, n-dimensional Riemannian manifold, and \et j : M ^^ N be 
the inclusion of a submanifold M, endowed with the metric g = j*7 inherited from 
N . If M is compact, then it is possible to construct a tubular neighborhood U 
around M; this is diffeomorphic to an open neighborhood of the zero section of the 
normal bundle of M, so there is a normal projection map a: U ^ M. In particular, 
there exists some ^0 > such that the set Msg , consisting of points in N whose 
Riemannian distance to M is less than 5o, is contained in U. (For details, see, e.g., 
Abraham and Marsden [1], Lang [23], Lee [25 •) Now, let jh : Mh ^-> iV be a family 
of inclusions of ?Ti-dimensional submanifolds M/j, parametrized by ft., each endowed 
with the Riemannian metric g^ = j;*7. If M^ lies inside the tubular neighborhood 
U and is transverse to a (i.e., Mh corresponds to a section of a), then it is possible 
to define the diffeomorphism iph = o\Mh '■ ^h — > M. 

An important case is when N — M", where n — m + 1 and 7 is the standard 
Euclidean metric, so that M C M" is an oriented Euclidean hypersurface. It is 
possible to define a signed distance function 6: C/ -> M on the tubular neighborhood, 
so that \S{x)\ = dist(a;,M) and \/S(x) — i^(x) is the outward-facing unit normal 
to M at a{x). Every point a; G C/ in the tubular neighborhood has a unique 
decomposition 

X ~ a{x) + 6{x)iy{x), 
so the normal projection map a: U ^f M can be written as 

a{x) — X — 5{x)v{x). 

Therefore, 

Va = I -V5®v- 5Vv = I -v®v- 5Vv = P + 6S, 

where P = I — iy(E)i^ is the projection map onto TM and S = —^v = — V^(5 
is the shape operator, or Weingarten map. (Note that Dziuk [TB], Demlow and 
Dziuk [13], Demlow [14j define a Weingarten map H = —S using the opposite sign 
convention, but this is less common in the differential geometry literature.) 

Instead of directly computing the tangent map Ta: U — >■ M, we can look at 
its adjoint, which "lifts" vectors on M to those on U. Given the pullback map 
a* : i}^{M) — > il^{U), the metric can then be used to identify covectors with vectors. 
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thereby obtaining a pullback map of vector fields X{M) — > X{U). Specifically, let 
Y e TyM and x e a^^{y) C U. Then define the lifted vector a*Y £ T^U satisfying 

X ■ a*Y = Ta{X) ■ Y, VX G T^U. 

In terms of the Riemannian sharp and flat maps, this can be written as 

[a*{Y)f' = a*{Y^) <^=^ a*Y = [a* {Y^)f . 

The following theorem allows us to compute this lifted vector explicitly, in terms of 
the signed distance function and shape operator. 

Theorem 4.3. Let M be an oriented, compact, m-dimensional hypersurface of 
]^m,+i ^j^^ ^ tubular neighborhood U. IfYG TyM and x S a~^(y) C U , then the 
lifted vector a*Y G T^U satisfies 

a*Y = {I + 5S)Y. 

Proof. Extend F to a constant vector field on K'"+^, so that Y — Vip{y) for the 
scalar function i^{x) = x ■ Y. Using the definition of the gradient V^" = (d'0) , 
and the fact that the exterior derivative d commutes with pullback, we have the 
following chain of equalities: 

a*Y = a* {W^p) = [a* (di/-)]* = [d (a* V)]" - V (V- o a) . 
Therefore, applying the chain rule, we get 

a*Y = Va(x) • V^- {a{x)) = {P + SS)Y={I + SS) Y, 
where the last equality follows from PY — Y. D 

Finally, when x G Mh, we can restrict to T^M^ by composing with the adjoint of 
jh, i.e., the projection Ph ^ I — Vh® Vh, which gives 

Yh=jla*Y = Ph{I + 5S)Y. 

This map j^a* = P^ (/ + 5S) is immediately seen to be the adjoint of the restricted 
tangent map Tiph = T (aJAf^) = T {a o jf^) = Ta o Tjh- In the next theorem, we 
bound the singular values of this map, thereby obtaining an estimate for the error 
term |j/ — J/i|| in the case of Euclidean hypersurfaces. 

Theorem 4.4. Given an oriented, compact, m-dimensional hypersurface M C M™"'""'^ 
with a tubular neighborhood U, let Mh be a family of hypersurfaces lying in U and 
transverse to its fibers, such that \\6\\^ jWi' ~ I'hWoo — > as /i -> 0. Then, for 
sufficiently small h, 



i/-Aii<c(ii5iL+ii^-^,iiy 



Proof. To obtain bounds on Y/j = Pha*Y, and hence on the singular values, suppose 
without loss of generality that |y| = 1. By the triangle inequality, 



|2 



|2 u^ |2 



l-\Yhr < l~\a*Y\' + |a*yr-|m 

For the first term, the eigenvalues of the shape operator are the principal curvatures 
Ki (x) , . . . , Kmix) for a surface parallel to M at x; as noted in Demlow and Dziuk 
|15) . Demlow [Tl], these are related to the principal curvatures at a{x) G M by 

Ki {a{x)) 



Ki{x) = 



1 + 5{x)ki {a{x)) 
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It follows that the eigenvalues of 5S at x can be estimated by 



d{x)Ki (a(x)) 



1 + S{x)K,i {a{x)) 
Since a*Y ^ {I + SS) r and |r 



1 + S{x)Ki {a{x)) 
1, this immediately implies 



<C\Six)\. 



l-\a*Y\ 



<C|<5| 



For the remaining term, observe that since Yh — Pha*Y, 



\Y,,\' = \a*Y - Vh iyn ■ a*Y)f = \a*Y\^ - {vn ■ a*Yf , 



and therefore 



\a*Y\^-\Y^,\^ ^ {i^h ■ a*Y) 



{Piy^ ■ a*Yf < \P,,h\' \a*Y\\ 



Now, 

\Pt^hf = W-^{v ■ i'h)\^ = l-{v ■ Vhf = (1 + i^ • i^/i) (1 - i^ • i^/i) < 2 (1 - i^ • Vh) , 



Vh\ , it follows that 

2 



and since 2{\ — v ■ Uh) = \ 

\a*Y\' - \Yu\'\ <\v- vnf \a*Yf <\v ^ v^f {\ + \\ - \a*Y\ 
Putting these together, we have 

\-\Yhf\<c{^b\ + \v^v^fy 
from which it follows that at each x G A/^, the singular values satisfy 
\\-a.,\<c{^8\^\v-v^f\, z = l 



<C\v-vA 



,m. 



Finally, applying (11), we obtain the uniform bound 



\S\\. 



\i^- Vh\ 



which completes the proof. 



D 



We now apply this theory to an important class of examples, where Mh corre- 
sponds to a family of piecewise- linear triangulations (as in Dziuk |16j , Demlow and 
Dziuk [I5j). or more generally, to the family of approximate surfaces obtained by 
degree-s Lagrange interpolation over each element of the triangulation (as in Demlow 
|14)). where the piecewise-linear case corresponds to s = 1. Here, the elements of 
this triangulation are assumed to be "shape-regular and quasi-uniform of diameter 
/i" [13] . Note that Mh is always constructed from an underlying piecewise-linear 
triangulation, even in the case of higher-order polynomial interpolation. Thus, by 
Corollary 4.2 we can define the V^ and Vr families of finite element differential 



forms on M/j, and obtain bounded cochain projections, even when s > 1. 

By Demlow Ti;, Proposition 2.3], for sufficiently small ft,, the surfaces Mh obtained 
by degree-s Lagrange interpolation satisfy 

(12) m^<Ch^+\ \W~Vh\\^<Ch\ 



Therefore, we obtain the following corollary to Theorem 4.4 



Corollary 4.5. If Mh is a family of surfaces approximating M, obtained by degree-s 
Lagrange interpolation, then \\I — Jh\\ < Ch'^^^. 
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Proof. Applying (12), we have 

11/ - All < C {\\S\\^ + \W- i^hWl) < Ch'+' + Ch"' < Ch'+\ 

which completes the proof. D 

This result generalizes Demlow [lU Proposition 4.1] — which applies only to scalar 
functions (k = 0) on hypersurfaces of dimension m = 2, 3 — to hold for arbitrary 
fc-forms, A: = 0, . . . , TO, on hypersurfaces of any dimension. In particular, the special 
case fc = 0, m = 2, s = 1, gives ||/ — Jft,|| < C/i^, which recovers the original estimate 
of Dziuk [TS] for piecewise-linear triangulation of surfaces in M.^. The correspondence 
between this framework, and that of Dziuk and Demlow, will be made explicit in 
the following worked example. 

Example 4.6 (The Hodge-Laplace operator on a 2-D surface). Let M be a closed, 
two-dimensional surface, embedded in M.^, and suppose the approximate surface Mh 
is obtained from degree-s Lagrange interpolation over a piecewise-linear triangulation 
Th- Assume that Th is contained in a tubular neighborhood of M, that its vertices 
lie on M , and that its triangles are shape-regular and quasi-uniform of diameter h. 

Take the continuous Hilbert complex to be the L^-de Rham complex on Af , i.e., 
W = L^^{M) and V — H^{M). Since Th is piecewise-linear and simplicial, we 
can take the discrete complex to be any of those considered in Arnold, Falk, and 
Winther [21 13]. For this example, let us take V^ to be the space of Vr /c- forms, and 
V,j ~^ to be the space of Vr+i {k — l)-forms. We emphasize that the fact that Th is 
a surface embedded in 'M? , rather than a flat region in M^, does not introduce any 
additional complications as far as the discrete complex is concerned. Indeed, the 
shape functions are defined with respect to a two-dimensional reference triangle, 
and this reference triangle can be mapped onto a triangle embedded in M? just 
as easily as one in M? . These shape functions can, likewise, be lifted up from Th 
to the curved triangles on the interpolated surface Mh- For nodal Lagrange finite 
elements (fc = 0), this observation was made by Dziuk [1,6] in the piecewise linear 
case, leading to the development of surface finite elements, while Demlow [H] later 
extended this argument to higher-order Lagrange polynomials. (Similar ideas had 
also been used in the development of isoparametric finite elements for Euclidean 
domains with curved boundaries.) 

Now, given some / G L'^Vt^{M), we obtain a solution {a,u,p) to the variational 
problem (Is]) on M. For the discrete variational problem (Is]), we can use the tubular 
neighborhood projection to take fh = oIm^/, thus obtaining a discrete solution 
{ah,Uh,Ph) on Mh- The modified discrete solution (crjj, ujj,^'^) — which is used only 
for the analysis, but is not necessary for computation — also lives on Mh, while its 
image {ih<7'i^,ihu'f^,ihP'h) hves on M itself. 

Therefore, assuming sufficient elliptic regularity, the "improved estimates" of 
Arnold, Falk, and Winther [51 13] yield the L? estimates for the modified problem, 

\\u-ihu'^\\^\\p-ihP'h\\<Ch-^'\\f\\H,.^,, 

||d (w - z.OII + ||a - *;,a;j| < C/i'- ||/IIh.-i , 

\\A{a-ih<j'h)\\<Ch^-^\\f\\^.^,, 
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which can be combined into the single estimate 



\W-ihUh\\ + \\p-ihPh\\ +h{\\d{u-ihu',^)\\ + \\(J - ih<Jh\\) 



Applying [Corollary 4.5 to account for the surface approximation error, we obtain 



the final error estimate for the discrete problem, 

\\u- ihUhW + \\p~ ihPhW +h{\\d{u~ ihUh)\\ + \\<J - ih<Jh\\) 

+ h'\\d{a~ z,c7,0|| < C (h^+^ WIWh^-. + h'^^ 11/11) . 

In particular, this implies that choosing isoparametric elements, with r — s, yields 
the optimal rate of convergence. 

The case k = Q and r — s — 1 corresponds to the lowest-order approximation 
of the Laplace-Beltrami equation, where Mh is piecewise-linear and V^ consists of 
piecewise-linear "hat functions" on Mh- In this case, the estimate above becomes 

\\u - ihUhW + lb - ^hPhW + h II V {u - ihUh)\\ < Ch^ ||/|| , 

which precisely recovers the estimate of Dziuk [TS], Demlow and Dziuk [15 . More 
generally, taking fc = with arbitrary r and s, we have 

\\u - ihUhW + \\P~ ihPhW +h\\^{u~ ihun)\\ < C {h-+^ ||/|l^,._i + h'+^ 11/11) , 

which agrees with Demlow [H] . 

On the other hand, we can also extend these estimates to the cases fc = 1, which 
corresponds to the mixed formulation of the vector Laplacian, and k = 2, which 
corresponds to the mixed formulation of the scalar Laplacian. For fc = 1, the 
estimate for general r and s becomes 

\\u- ihUhW + \\p- ihPhW +/i(||V X {u~ihUh)\\ + \\cr- ihCrhW) 

+ h^ ||V {a - z,f7,0|| < C (r+i WfWHr-. + h'+' 11/11) , 
while for k — 2, we obtain 

\\u - ihUhW + \\p - ihPhW +h\\a - ihO-hW 

+ h^\\V-ia~ z,f7,0|| < C (h^+^ Wn^r-. + h^+' 11/11) . 

4.4. Other variational crimes. The variational crimes framework developed in 
|Section"3| is quite general, representing a natural extension of the Strang lemmas 
from Hilbert spaces to Hilbert complexes. As such, the standard "crimes" that are 
typically analyzed in Hilbert spaces with the Strang lemmas may now be analyzed 
in the setting of Hilbert complexes. These crimes — including numerical quadrature, 
approximate coefficients, approximate boundary data, approximate domains, as 
well as isoparametric and other geometric approximations to more complex domain 
shapes — can all be represented as an approximate bilinear form Bh, an approximate 
linear functional Fh, and an approximation subspace Vh ^ T^, as in (|2|. In addition, 
techniques such as mass-lumping, which yield a number of benefits, such as discrete 
maximum principles and more efficient evolution algorithms for parabolic equations, 
are often analyzed in a similar way, and as such may now be analyzed in Hilbert 
complexes through the framework developed in |Section 3[ 
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5. Conclusion 



We began the article in [Section 2| with a review of the mathematical concepts 
that play a fundamental role in finite element exterior calculus, as developed in 
Arnold, Falk, and Winther [3^; these included abstract Hilbert complexes and their 
morphisms, domain complexes, Hodge decomposition, the Poincare inequality, the 
Hodge Laplacian, mixed variational problems, and approximation using Hilbert 
subcomplexes. In [Section 3[ we then considered approximation of a Hilbert complex 
by a second complex, related to the first complex through an injective morphism 
rather than through subcomplex inclusion. We developed several key results for this 
pair of complexes and the maps between them, and then derived error estimates 
for generalized Galerkin-type approximations of solutions to variational problems 
using the approximating complex; these estimates can be viewed as generalizing 
the results of Arnold, Falk, and Winther [3] to "external" approximations. Our 
main abstract results are thus essentially Strang-type lemmas for approximating 
variational problems in Hilbert complexes. 

As an application of the new framework developed in [Section "Bj wc developed a 
second distinct set of results in [Section 4[ for the case of the Hodge-de Rham complex 
of differential forms on a compact, oriented Riemannian manifold. We first reviewed 
Hodge-de Rham theory, and then considered a pair of Riemannian manifolds related 
by diffeomorphisms. We then established estimates for the maps needed to apply the 



generalized Hilbert complex approximation framework from Section 3| subsequently 
specializing this analysis to the case of a Euclidean hypersurface, with approximating 
hypersurfaces living in a tubular neighborhood. The surface finite element methods, 
as analyzed in Dziuk [T^, Demlow and Dziuk [TS], Demlow [T3], fit precisely into 
this class of approximation problems; as such, we illustrated how our results recover 
the analysis framework and a priori estimates of Dziuk J16j . Demlow and Dziuk 
[15], Demlow [14], and also extend their results from scalar functions on 2- and 
3-surfaces to general k-iorms on arbitrary dimensional hypersurfaces. Our results 
also generalize those earlier estimates from nodal finite element methods for the 
Laplace-Beltrami operator to mixed finite element methods for the Hodge Laplacian. 
By analyzing surface finite element methods using a combination of general tools 
from differential geometry and functional analysis, we are led to a more geometric 
analysis of surface finite element methods, whereby the main results become more 
transparent. 

There remain a number of interesting and challenging problems that were not 
addressed in the current article. One such problem is the extension of the pointwise 
error estimates of Demlow [14 for 0-forms to general fc-forms; this analysis relies 
on known results for the Green's function of the Laplace-Beltrami operator on 
the continuous surface (cf. [1]), and analogous results would be needed for general 
fc-forms. A second problem of interest is an extension of the Hilbert complex 
framework to more general Banach complexes, as would be needed to handle some 
nonlinear problems. This leads to a third interesting problem, which would involve 
the extension of the weak-* convergence and contraction frameworks, used for 
adaptive finite element methods for linear [S] [5S] and nonlinear [22! problems, to 
the setting of finite element exterior calculus, as well as to the surface finite element 
setting. 
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